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\^ \ Abstract 

^^ ' Wc consider a stochastic directed graph on the integers whereby a directed edge 

between i and a larger integer j exists with probabihty Pj-i depending solely on the 

Or . I distance between the two integers. Under broad conditions, wc identify a regenerative 

P I ■ structure that enables us to prove limit theorems for the maximal path length in a long 

(-H \ chunk of the graph. The model is an extension of a special case of graphs studied in 

[TB] . We then consider a similar type of graph but on the 'slab' Z x /, where / is a 

finite partially ordered set. We extend the techniques introduced in the in the first part 

of the paper to obtain a central limit theorem for the longest path. When / is linearly 

ordered, the limiting distribution can be seen to be that of the largest eigenvalue of a 

|/| X |/| random matrix in the Gaussian unitary ensemble (GUE). 

> 

^ ' 1 Introduction 

ly-^ . Consider a random directed graph with vertex ^ = Z, the integers. A pair of integers 

^D \ {hj) is declared to be an edge, directed from i to j, with probability Pj-i which depends 

only on the difference j — i, and this is done independently from pair to pair. We assume 

that pfc = for all fc < 0, so there are no directed edges from a larger integer to a smaller 

one. We are interested in limit theorems (law of large number and central limit theorem) 

^\ ' for the maximum length T[l,?i] of all paths from 1 to n, as n — t- oo. The problem as such 

j^ ■ is related to last-passage percolation. 

Unlike nearest-neighbour graphs |28^ [3]. the quantity T[l,n] does not have a direct sub- 
additive property. It turns out that, a related quantity, namely the maximum -L[l,n] of all 
paths in the restriction of the graph on {1, . . . , n}, has an almost sub-additive property (see 
([2|)) and thus L[l,n]/n — >• C, almost surely, for some deterministic constant C < 1. It is 
later shown that any two vertices are almost surely eventually connected by a path, and 
thus T[l,n] has the same asymptotic properties as L[l,n]. The minimal condition we need 
to carry out our programme is 



^(l-pi)---(l -pk) < oo. 



fc=i 
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Under this condition, we can identify a random subset .5^ (we call it "skeleton") of TL 
whose points form a stationary renewal process (see Sections [3] and |3|) over which the graph 
regenerates and has the property that any element v of =5^ is connected by a path (directed 
either towards v or away from it) to any other vertex in Z. The quantity -L[l,n] becomes 
additive over the regenerative set 5^ enabling us to prove, under the stronger condition 

oo 

^A;(l -pi)---(l-pfc) < oo, 
fc=i 

a (functional) central limit theorem. The latter condition implies finiteness of variance of the 
longest path between two successive points of .5^ . To prove the latter assertion, we provide 
a rather non-trivial algorithmic construction of the last non-positive element of .5^ . This 
construction is related to the so-called coupling-from-the past method for perfect simulation 
[331 fT9] and is the topic of Section [6] which is based on the properties of two stopping times 
studied in Section [5l The central limit theorem is proved in Section [3 

We then consider an extension of the random graph on the vertex set Z x /, where / is 
a partially ordered set under some partial order -< possessing a minimum and a maximum 
element. We let an edge from {x^i) to (y,j) exist with probability that depends on y — x and 
on i and j, and only when y — x > and i ^ j. We let L^v be the length of the longest path 
in the restriction of the graph on {0, . . . , A^} x i and show that the law of Lat, appropriately 
normalized, satisfies a functional central limit theorem such that the limit process {Zt, t > 0) 
is 1/2-self-similar, non-Gaussian, continuous process with Zi having the law of the largest 
eigenvalue of an a |/| x |/| random matrix in the Gaussian Unitary Ensemble (GUE) [2]. 

The case where all the pk are equal to p corresponds to a directed version of the classical 
Erdos-Renyi graph [3]. Indeed, let Gn,p be the Erdos-Renyi graph on the set of vertices 
{1, . . . ,n}. To each {i,j} which is an edge in Gn,p we give an orientation from i A j to 
i V j. The directed graph thus obtained is precisely the restriction of our graph on the set 
{1, . . . ,n}. This model was also studied in |18j . In this paper, we obtained, among other 
things, sharp estimates for the C = C{p) as a function of p. Besides purely mathematical 
interest, this model is motivated by applications in Mathematical Biology (community food 
webs) [311 El [30], in Computer Science (parallel processing systems) [22], and in Physics. 
Allowing the connectivity probability to depend on the distance between two vertices i and 
j means larger modelling flexibility on one hand while making the model more realistic on 
the other. 

In |18j we developed a generalisation of Borovkov's theory of renovating events [9l [TOl [TTl 
[T2I [13] in order to construct a Markov chain in infinite dimensions describing the "weights" 
of vertices. As a matter of fact, in [18], the random graph was a special case of a more 
general dynamical system (the "infinite bin model") with stationary and ergodic input. 
In this paper, we follow a different approach, one that is applicable specifically for cases 
where there is independence between links. In such a case, the approach has the advantage 
that it is more elementary using, essentially, renewal theory and coupling between renewal 
processes. 

2 The line model 

We are given a set of numbers {pj , j € N) , such that 

0<Pj<l, j G N. 



and consider {oij, i,j £ 7j, i < j) as a collection of i.i.d. random variables with common 
law 

P(ao,i = 1) = 1 - P(ao,i = -oo) = pj-i. 

Based on this collection, we build a directed random graph G on Z with edges 

E = {{i,j) eZxZ: i<j, aij = 1}. 

We shall occasionally refer to the restriction G[i,j] of the graph on the vertex set {i,i + 
1, . . . , j} (deleting all edges with either of the endpoints not in this set). We are interested 
in the behaviour of longest paths. A path vr is an increasing sequence of vertices vr = 
{io,ii, ■ ■ ■ ,ii) successively connected by edges, i.e. ajo,ii = • • • = ai^_^ j^ = 1. The number 
^ = |7r| of edges is the length of this path. 

For any i > 1 and any increasing sequence (io,ii, . . . ,i^) of vertices we conveniently 
define 

\{io,ii,...,ii)\ = (aio,Ji +aji,J2 H \- ct^.t,ie)~^ ■ (1) 

Clearly, this quantity is if one of the summands takes value — oo; otherwise, it equals i. 
In other words, |(zo, ii, . . . , i^)! > if and only if (io, h, ... ,11) is a path. 

We say that there is a path from i to j ii Iq = i, i^ = j; we denote this event hy i -^ j 
and may also express it by saying that i leads to j or that j is reachable from i. 

We let T[i,j] be the maximum length of all paths from i to j. Unlike nearest-neighbour 
directed graph models (see, e.g. [2^), this quantity does not have a subadditivity property. 
To remedy this we let L[i,j] be the maximum length of all paths from some i' > i to some 
/ < j, i-e., 

L[hj] = . max T[^''/]- 

That is, L[i,j] is the longest path of the restricted graph G[i,j]. Clearly, L[i,j] has the 
same law as L[0,j — i]. It is also clear that L[i,j] is subadditive in the sense that 

L[i,k]<L[i,j]+L[j,k] + l, i<j<k. (2) 

Indeed, if vr is a path of maximal length in G[i, k] then its restriction vr' on G[i, j] has length 
at most L[i,j] and its restriction vr" on G[j,k] has length at most L[j,k]. Now the length 
of TT is equal to the length of vr' plus the length of vr" plus, possibly, 1, if j is not a vertex 
of vr. By the subadditive ergodic theorem [25l p. 192], there exists a deterministic C £ [0, 1] 
such that 

F{limL[i,j]/j = G) = l. (3) 

Some of the results below do not depend on the independence assumptions between the 
random variables aij . It is often necessary to define the model on an appropriate probability 
space. We do this as follows. Let S = {Sj,j G Z) be a collection of independent {—00, 1}- 
valued random variables with 




F{6, = 1) = 

Let d^^' , i S Z be i.i.d. copies of S. The probability space Q consists of u = (5*-*^ , i G 
The random variables Oij are then defined by 

aij(w) = 6\%. 



The sigma-field is the standard product signia-field. A natural shift 9 on il is the map 
defined by 

a; = (i^^«)^^w=(i^5(^+^)). (4) 

Hence 

aij{9uj) = 5^-^P = <5[]+;5-(i+i) = ai+ij+i(^)- 

The random variables L[i,j] are all defined explicitly on VL via L[i,j] = maxj<j,j<j^<...<j^<j |(io, ■ ■ ■ ,11) 
where (iq, . . . , ii) is the random variable defined by ([T]). It is in this sense that the law P of 
the model is ^-invariant on Q. Moreover, 9 is ergodic. In fact, the result that the asymp- 
totic limit of L[l, n]/n exists depends only on this 0-invariance, so it holds for more general 
models where the law of d is not that of independent random variables. 

A word on notation: If {An,n G Z) is a collection of events of ft and r is Z- valued random 
variable on Q then A^ denotes the event containing all w G il such that lo E ^r{w)- 

3 The skeleton 

For the purposes of this section, let be the space defined above, 9 the natural shift (j3]), 
and let P be a ^-invariant probability measure. In addition, assume that 9 is ergodic, i.e. 
that the invariant sigma-field is trivial. Recall the shorthand {i -^ j} = {T[i,j] > 0} for 
the event that there is a path from i to j. Consider, for each n £ I^, the events 

A'^ '■= \\{n -^ j} = {any j > n is reachable from ?i} 

j>n 

A~ := (I {j ~-> 7i} = {n is reachable from any j > 7i}. 
j<n 

The following is an immediate consequence of the definitions: 

Lemma 1. (i) The sequence ((j4~, A^),n G Z) is stationary and ergodic. (ii) For each n, 
the events A~ and A^ are independent and F{A^) = F{A~) = F{Aq). 

We are interested in the random set 

y{uj) :={n£Z: u £ A+ n A'}, (5) 

and refer to it as the skeleton of the random graph. The terminology is supposed to be 
reminiscent of a point of view described next. 

Let ^(E) C ZxZbeapartial order (i.e. if (i,j),(i,^) e ^(£^) then (i. A;) G ^(£^)) which 
contains the set of edges E. In fact, take 0^{E) to be the smallest such set. Necessarily, 
^{E) = {(i, j) S Z X Z : i -^ j}. A subset [/ of Z is totally ordered under the partial order 
-^ if for any distinct i,j £ U we either have i -^ j ox j -^ i. We say that a totally ordered 
subset U is special if it has the stronger property that for all distinct i, j with i £ U and 
j £ V, we either have i -^ j 01 j -^ i. Clearly, the union of special totally ordered subsets 
is special; thus we can speak of the maximal special totally ordered subset and we refer to 
it as the skeleton of the partial order. Adopting this definition, it is now clear that the set 
y defined by ([5]) is the skeleton of the partial order -^ on Z. In [1] the elements of 5^ are 
referred to as posts. In fact, [T] uses J^ in order to derive limit theorems of the number A'^ 
of linear extensions of the the random partial order ^^ on {1, . . . , n}. 

For a general partially ordered set, a skeleton may not exist. However, in our case, the 
condition P(A^ n A^) > is sufficient for ^ to be almost surely infinite. 



Lemma 2. If X := ^{Aq n Aq—) > then .5^ is an a.s. infinite set. 

Proof. Let 6 be the shift defined by dH). Then, for all w, y{uj) = y{9uj). Since P is 
^-invariant, the result follows. D 

Assuming that A = IP(^o^ H Aq) > 0, we may then, equivalently, consider ^S^ as a 
stationary-ergodic point process on the integers with rate A because A = P(0 G o?*'). We let 
r„, n E Z be an enumeration of the elements of ^ according to the following convention: 

• • • < r„i < To < < Ti < La < • • • 

In particular, Fq is the largest non-negative element of J/'. 
We can now strengthen the subadditivity property ([2]) for L: 

Lemma 3. For all integers m < n, 

Proof. To see this, consider the interval [ri,r„] and a path vr* of length L[ri,r„]. Then 
this path must visit all the intermediate skeleton points Ti, . . . ,r.„. Lideed, suppose this 
is not the case and vr* does not visit, say, F/, for some 1 < I < n. Consider an edge {i,j) 
belonging to vr*, with i < F; < j. By the definition of F^, both (i, Ti) and (F;, j) are edges of 
the random graph G. Therefore we can increase the length of vr* by 1 if we replace the edge 
{i,j) by two edges {i,Ti) and {Ti,j). This leads to the contradiction since vr* has length 
L[Fi,F„] which is, by definition, maximal. D 

4 Regenerative structure 

Throughout, we make use of the following two conditions: 
[CI] < pi < 1 

oo 

[C2] ^(l-pi)...(l-p,)<oo. 
fc=i 

We also sometimes write Qj = 1 — pj. For each j G Z we consider its immediate neighbours: 

ri{j) := min{A; > j : aj^k = 1} 

?(j) := max{i < j : atj = 1}. (6) 

See Figure 1. The distances of these vertices from j are denoted as follows: 

vU) ■■= vij) - j 

m-=j-m- 

Notice that {C{j),j G Z) and {rj{j),j G Z) are identically distributed sequences, and that 
each one is a sequence of i.i.d. random variables. Furthermore, for each j G Z, 

iCiJ + 1), e(i + 2),...)AL {7]{j - 1), 7?(i - 2), . . .) 



[W)) ooooMooo riQ) 



Figure 1: Notation used: £,{j) is the first vertex below j that is connected 
to j; correspondingly, rj{j) is the first vertex above j connected to j. 

Henceforth, we shall let ^ be a random variable with distribution the common distribution 
of C(j) and r]{j): 



>n)= P(^(0) > n) = P(?7(0) > n) = (1 - pi) • • • (1 - p„,), n G N. 
Define next the events 

V v—l 

j=u+l j=u 



4+ ~i 4+ A~ -^ A~ 

lim Kv = ^t, l™ Kv = A-. (8) 



for which, clearly, 
with 

Furthermore, 

Kv n Kw ^ Kw^ ifu<v<w, (9) 

a property we shall use in Section [H Observe also the following: 
Lemma 4. For all integers u < v, 

V j — 1 V 

j=u+l i=u j=u+l 

v—l V v—l 

j=ui=j+l j=u 

j>u i=u j>v 

V 

j<v i=j+l j<v 

Proof. We prove the first equality. That A+„ C ^l^^^_^_l U^~^ {i -^ j} is immediate from the 
definition ([7]). To prove the opposite inclusion, assume that u > v + 1 (otherwise there is 
nothing to prove) and that for all integers j € [u+l, f ] there exists an integer i G [u, j — 1] such 
that i -^ j. Fix j > u and pick ii to be the largest among the vertices between u and j — 1 
such that ii -^ j; necessarily, Oj^j = 1. Then pick the largest vertex Z2 among the vertices 
between u and ii — 1 such that «2 '^ ii, and continue this way. Since ii > i2 > ■ ■ ■ > u, 
it follows that this process terminates with some ik = u. Since {u = ik,ik+i, ■ ■ ■ ,ii,j) is a 

6 



path, we have that u -^ j. The second equality for A^^ now follows from the definition ^. 
The relations for A~^ follow similarly. The third (respectively, fourth) line is obtained by 
sending v to +00 (respectively, u to —00) in the first (respectively, second) one. D 

This lemma tells us that A^^ is the intersection of -y — n independent events. Indeed, 
since ^(j) = j — ^{j) we have 

Kv = {^(^ + 1) < 1, ^(u + 2) < 2, . . . , e(f ) <v-u}, (10) 

and the random variables ^(n + 1), . . . ,£,{v) are i.i.d. Similarly, for A~^, 

A-^, = {7?(u) <v-u,...,ij{v-2)<2,r]{v-l)< 1}. (11) 

Moreover, since 

(e(n + 1), e(n + 2), . . . , Civ)) ^ iviv - 1), 7?(r; - 2), ... , r]{u)) 

we have that ¥{A^^) = ¥{A~^). Similarly, both A^ and A~ are intersections of infinitely 
many independent events: 



A+=f]{C{j)<j-n} (12) 

j>n 



j<n 

and P(j4+) = F{A^). The skeleton 1^ can be expressed as follows: 



y = {neZ: supfj{i) < n < inf ^{j)}. (14) 



Regarding ^5^ as a point process, we see that it has rate 

00 00 

A = p(o e ^) = nKf = iU^im < j))' = nti - ip(^(o) > ^■)]'- 

Since 

P(e(0) > j) = P(ao,i = • • • = aoj = 0) = (1 - pi) • • • (1 - p^), (15) 

we have 

A=n[i-(i-pi)---(i-rf (16) 

and so 

[C2] ^^ A > ^^ E[^(0)] < 00. 

Consider now two successive skeleton points T^ and r^+i and let '^^(a;) be the restriction 
of w on [rfc,rfc+i): 

^fc:=(d("), rfc<n<rfc+i), kGZ; 

we refer to it as the fe-th "cycle" . We next show that the sequence of cycles have a regen- 
erative structure in the following sense: 

Lemma 5. The cycles C^k^k G Z) are independent and ('^fc,A; G Z — {0}) are identically 
distributed. In particular, the skeleton vertices (F^, A; € Z) form a stationary renewal process. 



Intuitively, Lemma [5] is based on the following observation. Suppose that is a skeleton 
vertex (i.e. condition on the event Aq n Aq). Then ^(1) > 0, ^(2) > 0, etc. In other words, 
^(1) = 0, ^(2) G {1, 2}, f (3) G {0, 1, 2}, etc. To determine the location of the next skeleton 
vertex after we need to find the first vertex j > such which is connected with every 
vertex between and j — 1. This means that, conditional on being a skeleton vertex, the 
location of the first skeleton vertex larger than does not depend on the {S^^\n < 0). 
Proof of Lemma O 

For k > 1, let ^^ be the sigma-algebra generated by {S^^\ . . . ,d^'''>) and let ^j7 be the 
sigma-algebra generated by {S^~^' , . . . ,S^~^'). It suffices to prove that, for any k > — 1, 
/ > 1, and any Bj: £ ^^", B^ £ ^,+ , 

p(r_i = k,B^,ri = i,B+ I To = 0) = P(r_i = k,B^ \ Tq = o)P(ri = i,b+ | Tq = o). 

Assume that Fq = (i.e. is a skeleton vertex). Then, by (J14p . 

..., fj{-2), 7?(-l)<0<ai), ?(2), ... 
In view of the latter inequality, we have 

r_i = max{n < : l^-nA+ = 1} 

= max{n < : . . . , ri{n - 2), fj{n - 1) < n < ^(?i + 1), ^(?i + 2), . . .} 

= max{n < : . . . , r?(n - 2), r/(n - 1) < n < ^(n + 1), ^(n + 2), . . . , ^(0)} =: f _i, 

where the last serves as a definition of a new random variable r_i. This random variable 
is ^"-measurable, where ,^~ is the sigma-algebra generated by {S^^\k < 0). Similarly, we 
define 

fi := min{n > : r?(0), . . . ,f?(n - 1) < n < f(n + l),^(n + 2), . . .}, 

a random variable which is ^"^"-measurable, where ,^^ is the sigma-algebra generated by 
{S^'^' , k > 0), and observe that, on {Tq = 0}, the random variables Fi and Fi coincide. Note 
that ^~ and ^"^ are independent. Hence, for k < —1, i > 1, we have 

P(F_i = A:,B-,Fi = /,B+|Fo = 0) 

_ P(F_i = k, S-, Fi = I, B+,A+ n Aq) 



l+nylo") 
p({f _i = k}nA^nBj:, {f i = ;} n ^+ n b+) 

P(A+)P(Ao-) 

P({f_i = fc} n ^^ n Bj;) P({fi = l}nA+n B+) 

nAtMAo) 
' -i=k,B- I A~)F{fi=l,B+ \A+). 



Note that 



P(F_i = k, B^ I Aq) = P(F_i = k, B~ I Aq n A+) 

= P(F_i = k, B~ I A~ n A+) = P(F_i = k, B- |Fo = 0). 



Similarly, 



P{Fi = 1,B+ I A+} = P{Fi = 1,5+ |Fo = 0}. 

D 



Corollary 1. The bivariate random variables 

(ri-Fo, L[ro,ri]), {T2-Ti, l[t,,T2]), ... 

are i.i.d. 

5 Two stopping times 

In this section, we study properties of the fohowing two random variables: 

/i := mi{i > : 1 ,- =0} 

u := infji > : 1.+ = 1|. 

These random variables are important in the algorithmic construction of Section [6l 

Note that — z/ is the first vertex < with the property that every vertex in the interval 
(— z^, 0] is reachable from —i/: 

V = inf{i >0: — z^-^0, — i/~^ —1, . . . , — i/ -w —i/ + 1}. 

Also, — /i is the first vertex < such that is not reachable from — //: 

fi = inf{i > : - i -/^ 0}. 

We will show that ^ is a defective random variable, i.e. that P(/x = oo) > 0, with 
conditional tail P(yu > n|/i < oo) comparable to the integrated tail of ^. We will also show 
that v is an a.s. finite random variable with the same number of moments as S,- 

Note first that both n and v are stopping times with respect to the filtration (^^, A; < 0). 
Observe that 

{// = oo} = n^:^,o = ^^- (17) 

Since condition [C2] is equivalent to F{Aq) > 0, we have 

¥{H = oo) > 0. 
On the other hand, 

oo 

{z. = oo} = n {A\^,r, 

n=l 

and, as we shall see below, this event has probability zero: 

P(z. = oo) = 0. (18) 

Let us first focus on the law of fi, conditional on {fi < oo}. This can be computed easily, 
from the definition of fi, and equations ([TT]) . (fT7|) . and ([15]). 

P(n < (1 < oo) = P(77(— A;) < k for all 1 < A; < n) P(r/(— tti) > m for some m > n) (19) 

n / oo \ 

= n ^P'(^(-^) < ^) ( 1 - n ^ivi-m) <m)\ 

fc=l ^ m=n+l ^ 

= (1 -9i)(i -9i92)---(i-gig2---9n)( 1 - n (1 -^i^2---gm) j 

^ m=n+l ^ 

Conditional on {fi < oo}, the random variable /x has a tail comparable to the integrated tail 
of C: 



Lemma 6. Suppose that [CI] and [C2] hold. There exist constants < Ci < C2 < oo such 
that, for all n > 0, 

Ci^F{^ >m) <F{fi>n\ n<oo) <C2Y^ P(^ > m). 

m>n m>n 

Proof. Since pi < 1, we have A < 1 (see (fTB]) ) and so 

P(/x < 00) = 1 - A^/2 > Q^ 

Using (fT9]) we have 

_. 00 ^00 

P(/i > n I ^ < 00) < Y^^^ X] ^(^M > "i) = ^ _ ^1/2 I^ '^'(^ > "^)- 

m=n+l m=n+l 

Hence Ci = 1/(1 — A^'^). To obtain a bound from below note that the first term on the 
right of p9|) is > P(// = 00) and so 



P(n < /x < cx)) = A^/2 I 1 - PJ P(r/(-m) < n 

^ m=n+l 

> A^/2jl-exp(- ^ P(C>m))) 

^ m=n+l "^ 

00 

m=n+l 

where g{x) = (1 - e-^)/x. Hence C2 = g{E^)X^/^/{l - A^/^)). D 

We next prove something stronger than l\18\\ . namely that z/ has the same number of 
moments as ^. 

Lemma 7. IfK^^ < 00 for some r > 1 then Ku^' < 00. 

Proof. By the definition of u and equation (|1U|) we have 

z^ = inf{n> 1 : ^(0) < n,^{-l) < n - 1, . . . ,^{-{n - 1)) < 1} 

Define a sequence of non- negative random variables xq, xi, X2, • • • by xq = and 

Xn = maxU(O) - n, ^(-1) - (n - 1), . . . , ^(-(n - 1)) - 1}, n > 1. 

Then 

u = inf{?i > 1 : x.„ = 0}. 

The Xn satisfy 

Xn+i = m&x{xn,S,i-n)) - 1, n > 0, 

and, since the S,{—n) are i.i.d., {xn,n > 0) is a Markov chain in Z+. We now make two 
observations that imply the statement of the lemma. First, if rc„ > J^ > then 

xn+i -Xn = {i{-n) - x„)+ - 1 < (e(-n) - K)+ - 1. 
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But E[(^ — K)^] < 1 for sufficiently large K. Therefore, after the Markov chain leaves the 
interval [0, K] (for sufficiently large K) it is majorized from above by a random walk with 
increments distributed like {^ — K)~^ — 1 whose mean is negative. By standard properties 
of random walks this implies that the return time Tx to the set [0, K] satisfies ET^ < oo 
if E((^ — K)~^ — ly < oo; and the latter is equivalent to E,^^ < oo. The second observation 
is that the Markov chain (x^) returning to the set [0, K] eventually hits point after a 
geometric number of trials. D 

Corollary 2. // [C2] holds then Kv < oo. 

6 Algorithmic construction of Fq 

In this section we give a method for constructing a specific skeleton point, e.g., the 
first one which is to the left of the origin. This is the point Tq. Besides the theoretical 
interest, such a construction will be used later for proving a central limit theorem; it can 
also be used in connection to a perfect simulation algorithm for estimating the value of 
C = lim„_^oo L[l,n]/n (see remarks at the end of the section). 

The idea for the construction of Tq is this: recall that —i^ which is the first vertex < 
which is connected to every point between —v and 0. We check whether —i/ is also reachable 
from every point from the left. If it is, we declare that —i^ is a silver point and stop the 
procedure. If not, there is a first vertex before —v which fails to be connected to — z/. Using 
the shift operator 6 defined in (|3|, this vertex is at distance fj,o0~'^ from —i/; in other words, 
this distance is the functional fi applied to the shifted uj, when the origin is placed &t —v. 
We then set ^[1] = u + jio9^^ , which is the location of the previous vertex, and z^[l] = v and 
this finishes the first step of the procedure. 

The second step of the algorithm is similar to the first one: we search for the first vertex 
— 1^[2] before — /^[l] which is connected to every vertex between — z^[2] and — z^[l]. We know 
that we can find such a vertex with probability one. If it also happens that — z^[2] is reachable 
from any point from the left, we stop and declare — ^^[2] as our silver 'pova.i. Otherwise, there 
will be a first vertex, — /u[2] < — 1^[2] which fails to be connected to — 1^[2]. 

The procedure continues in the same way, until the first silver point is found, and it 
will be found with probability one. This first silver point will have the property that it is 
reachable from every point from the left and is connected to every point up until the origin; 
see Lemma [10] below. The distribution of this first silver point is well-understood and this 
is the content of Lemma [9l In fact, we will show that there are infinitely many silver points 
which form a (delayed) renewal process backwards; see Lemma [T2l Finally, in Theorem [1] 
we show that among the infinitude of silver points we can pick a gold one, namely the point 
To. 

To define the algorithm explicitly, we consider a sequence of N U {+oo}-valued stopping 
times relative to the filtration {^k,k > 1), defined as follows. Let 

u[l] := V 

^[1] := 1/ + ^o^-'^ = inf{j > 1/ : 1^- =0}, (20) 
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and, recursively, for k >2, 

u[k] := mf{j > fi[k - 1] : 1^+ = 1} 

/i[A:] :=z^[A;]+^o0-^W =inf{j>z^[A:] :1,- =0}, (21) 

where 9 is the natural shift Q. It is understood that if for some k we have ^i[k] = oo then 
i/[j] = /_i[j] = oo for all j > k + 1. We thus obtain an increasing sequence of stopping times 

i^ = zy[l] < ^[1] < z^[2] < /i[2] < z^[3] < ^[3] < • • • 

which (since P(/u = oo) > 0) is eventually equal to infinity. It is convenient to think of these 
stopping times as the points of an alternating point process (the /i-points and the z^-points). 
In words, the sequence of these stopping times is defined by first laying a i/-point in location 
z>[l]. Then, as long as r/(— (i/[l] + i)) < i for i = 1,2,..., we place no point in location 
v[l] + i. At the first instance i at which r/(— (i/[l] + i)) > i, we place a /.t-point in location 
v[l] + i and call it fi[l]. The random variables (77(— (z^[l] + i)),i > 1) are independent of 
z^[l], and so the event that we place a ^-point in a finite location is independent of z^[l] and 
has probability F{fi < oo) = 1 — A^'^. The procedure continues in the same way: having 
placed /^[/c] < oo, we decide, independently of the past (i.e. '^~ij.]) whether to create a new 
/x-point or not (i.e. place it at infinity). If we do create a new /x-point //[/c] then, clearly, 
z^[A: + 1] is also finite and z^[/i; + 1] — vlk] has the same distribution as z/[2] — ^[1] conditional 
on fi[l] < oo. Thus for each to, the recursion stops at the index 

K := infj/c > 1 : n[k] = oo}. (22) 

From the discussion above we immediately obtain: 

Lemma 8. Assume that [CI] and [C2] hold. Then K is a geometric random variable with 

¥{K >k) = {l- A^/2)fc^ ]^ > 0. 

By definition, ^[K] = oo but fi[K — 1] < oo. Hence 

I'lK] < oo, a.s. 

Note 1. We stop for a minute to point out that the whole purpose of the construction of 
these random variables is the random variable i'[K]. In other words, for each uj £ il., we 
apply recursion (I20p - (l2ip to obtain the alternating sequence of v and /i- points, through 
them we define that index K as in ()22p and, finally, vlK]. Thus, I'lK] is a well-defined 
(measurable) function of lo. We refer to —v[K] as the first silver point before 0. 

Although K depends on the whole alternating process (z^[A;], /-f[A:]), k > 1), we can identify 
the law of i^[K] as follows: 

Lemma 9. On a new probability space, let K,ijji,ijj2,ip3, ■ ■ ■ be independent random variables 
with distributions 

P{K>k) = {l-\^/^f, k>0 

ipi = {v[2\ - v[l] I ^[1] < oo) = ( inf{j > ^f : 1^+ = 1} | ^ < oo) , i > 2. 
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Then, assuming [CI] and [C2], 

K-l 

u[K] = iJi + ^ il^i+i. (23) 

Proof. It follows from 

K-l 

using a simple probabilistic argument as described above. D 

The reason we are interested in the random variable I'lK] is the following: 

Lemma 10. Assume [CI] and [C2] hold. Then for ¥-a.e. co 

^^^Z,lK]^At^[K],o- (24) 

Note that replacing the index n in a sequence of events An by a random index N amounts 
to defining the event A^ = {oj £ Q : there exists n such that n = N{lj) and co G An}. 

The meaning of (I24p is that the vertex I'lK] of the random graph has the property that 
there is a path from every j < iy[K] to i^lK] and there is a path from i^lK] to every i such 
that i^lK] < i < 0. Our goal is to identify a skeleton point. Whereas i^[K] is not a skeleton 
point for sure, there is a positive probability that it is. 

Proof of Lemma\l(A If K = k, for some k > 1, then fi[k] = oo but fi[k— 1] < oo, so vlk] < oo 
and 1^- = for all j > z>[A;]. Hence 



-j-i^m 



{K = k}c n ^:,,-.w = ^: 



j>u[k] 

by ([8|). Also, ii K = k, then i^[/i;], z>[A; — 1], . . . , v[l] < oo and so 

{K = k} c ^_j,[^,] _^[fc„i] n ^_^[fc_i] _^,[fc_2] n • • • n ^_^,[i]^o ^ "^-z.[fc],o' 

by ([9]). But iiT is a geometric random variable and hence K < oo, a.s. D 

We also have the following result concerning moments of I'lK]: 

Lemma 11. Assume [CI] and [C2] hold. If, in addition, there exists r > 1 such that 

EC^^ < oo, then Ei^[KY < oo. 

Proof. We have that Ei/[iir]'' < oo if Ez^'' < oo and E{fi^'\fi < oo) < oo. The latter holds if 
]g^r+i ^ ^^^ ^j^^ ^]^^g jg ^ simple consequence of Lemma [G] On the other hand, Ku^ < oo 
holds if E^'' < oo, as proved in Lemma [T] D 

Whereas [CI] and [C2] imply P(i/[ir] < oo), we need finite variance for (^ in order that 
we have finite expectation for I'lK]. 

We next construct a further sequence of stopping times. 

a[l] <a[2] <■■■ 



13 



as follows. Assume that [CI] and [C2] hold. Recall that the random variable i^[K] is a.s. 
finite; it maps O into N. Hence we can define u[K]o9"' for any n ^TL and also v\K\oQ-^ for 
any measurable J : $7 — t- Z. We define o"[j], j > 1, recursively: 

o{\\ = i.[K] 

c7[j + 1] = a[j] + i^[K]oe--^^\ j>l. (25) 

Intuitively, given u, we first construct i^[K] by (f20]l - (f2T]l and place a point a[l] at i^[K]. 
We then shift the origin to —u[K] and repeat the recursion with uj' = 6~'^^^'{uj) in place Q| 
of uj, thus obtaining a new random variable, u[K]o9~'^^^'. We place another point a[2\ at 
distance G i'[K]o6~'^i^^' from cr[l]. The procedure continues in the same way. We refer to 
— 0"[1], — o"[2], ... as the sequence of silver points. 

Lemma 12. Assume that [CI] and [C2] hold. Define the point process with points crlj], 
j > 1, as in (]25p . This is a renewal process on N, i.e. the random variables (t[1], a[2] — a[l], 
(7 [3] — a[2], . . . are i.i.d. with common distribution ([23]) . 

We are now ready to construct the first gold point Tq. 

Theorem 1. Assume that [CI] and [C2] hold. Define the sequence {i'[k],fi[k],k > 1) through 
(|20l) - (f2T]) which is used to define the random variable i^[K]. Based on this, define the sequence 
(o"[j],i > 1), through (|25p . In addition, let 

M := sup{^(z) - i}, 

J := inf{j > 1 : a[j] > M}. 

Then 

To = -a[J]. 

Before proving the theorem, let us observe that the random variables defined in the 
theorem statement are a.s. -finite. By [C2], i.e. that E^ < oo, implies M < oo, a.s. 

P(M > m) = P(^(i) -i>m, for some i > 1) 

oo 

< ^P(C(i) > i + m) (26) 

i=l 

CXD 

< Yl nm>i)<^i- (27) 

i=m+l 

By standard renewal theory, it is easy to see that J, the first exceedance of M by the random 
walk (cr[j],j > 1), is also a.s. -finite and hence (j[J] is an a.s. -finite random variable. 

Proof of Theorem [I] Owing to Lemma [JOj we have that 

for all j EN, a; e ^I^y, n ^+^^,^0, P-a.e. a; e Q. (28) 

Also, 

{M < a[J]] = {e(l) < a[J] + 1, i{2) < a[J] + 2, . . .}. (29) 
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Fix n G N and observe that, from the definition of M and the expressions (jlOp . ()12p for 
^i„ ^^^'^ ^-n' respectively, 

AZr, n Al„ n {A/ < n} = Ai„ n ^l„^o n {^(i) < i, e(2) < 2, • • • } 

= AZn n U(-" + 1) < 1, . . . ,e(0) < n,^(l) < 1,C(2) < 2, • • • } 

= A_^ n A_„ 

= {n G y}. 



Combining this with (128p and (I29p we obtain 

-o-[J] G =y, a.s. 

It is clear, from the algorithmic construction (I20p -()2ip of the sequence {v[k], ii[k\,k > 1), 
from the algorithmic construction (j25p of the (o"[j],j > 1), and the definition of J, that 
there can be no point of ^ between — cr[J] and 0. Therefore — o"[J] is the largest negative 
point of y. D 

Remark 1. Possible extensions: The algorithmic construction proposed above may be used 
in a general stationary ergodic framework. In particular, one can easily generalise first-order 
results (the functional strong law of large numbers). Under reasonable assumptions, one can 
again prove the finiteness of ^(0). This will imply the finiteness of 7?(0) and, in turn, the 
existence of the stationary skeleton. Then the functional strong law of large numbers will 
follow using well-known tools. 

Remark 2. Simulation and perfect (exact) simulation of the value of the limit C: This 
depends in a complex way on an infinite number of variables, and one cannot expect an 
analytic closed form expression. But one can estimate it by running a MCMC algorithm. 
One can also use the regenerative structure of the model to run the simulation in backward 
time using the idea of "cycle-truncation" that leads to a simple implementation scheme; 
c.f. [20] for more details However, each such an algorithm gives a biased estimator of the 
unknown parameter, in general. 

In [18], we considered the homogeneous case {pj = p, for all j). In particular, in [181 §10] 
(see also [181 §4] for theoretical background), we obtained a stronger result by proposing 
an algorithm for the perfect simulation of a random sample from an unknown distribution 
whose mean is the limit C under consideration. The standard MCMC scheme provides an 
unbiased estimator for this limit. 

The ideas behind that algorithm may be efficiently implemented in a number of similar 
models, e.g. in models with long memory ( see, for example, [IS]). In fact, in [18] . we 
developed the algorithm for a more general model (we called it "infinite-bin model") and 
under general stochastic ergodic assumptions. 

7 Central limit theorem for the maximum length 

Assume now that [CI] holds and 

[C3] ^fc(l-pi)...(l-p,)<oo. 

fc=i 
From (|15p we see that this is equivalent to 
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[C3'] E^^ < oo. 
Lemma 13. //[CI] and [C3] hold then K\Tq\ < oo. 

Proof. By Theorem [H |ro| = cr[J] = min{cr[j] : j > 1, a[j] > M}. Recall that a[l] < 
a[2] < • • • are points of a renewal process. This renewal process is clearly independent of 
M = supj>]^{^(i) — i}. By standard renewal theory, Ecr[J] < oo if EM < oo. But the tail 
of M was estimated in ()27p . The same inequalities now show that E^'^ < oo is sufficient for 
EM < oo. D 

The maximum length L„ of all paths from some i >0to some j <n satisfies the following 
central limit theorem. 

Theorem 2. Suppose [CI] and [C3] hold. Let 

a^:=Y&i{L{Ti,T2]-C{V2-Ti)). 

Define 

L\„i] — Cnt 
in{t) := "% ^ , t>0, nGN. 



Then the sequence of processes £n, in the Skorokhod space D[0, oo) equipped with the topology 
of uniform convergence on compacta Q/, converges weakly to a standard Brownian motion. 

Proof. By Lemma [TBI we have E|ro| < oo. Hence ELi < oo. But the r„ form a stationary 
renewal process. Therefore, ELi < oo implies that the variance of r2 — Li is finite. Since 
L(ri, r2] < r2 — Li, we have a^ < oo. The constant C, defined as the a. s. -limit of Ln/n-see 
([3]), is also finite and nonzero. Lemma[3] shows that (L„, n > 0) is a (stationary) regenerative 
process. The result then is then obtained by reducing it to Donsker's theorem. This is 
standard, but we sketch the reduction here for completeness. Let ^n be the cardinality of 
Jy^ n [0, n] (the number of Tj in the interval [0, n]): 

$„ := \yn[0,n]\ = ^1(0 < Tj < n). 

So r$^ < n < r$^+i. Write 

L[^t] = {L[nt] - Lr^^^^^} + Lr^^^^^ 

The quantities in brackets on both lines are tight and so they are negligible when divided 
by x/n. So instead of £n(t), we consider 

tn{t) := 'r,, ^'' = %r^ + -TTj^ j:{L(r._i,r,] - c(r, - r,_o} (so) 

The last term is the one responsible for the weak limit of in (and hence of £„). To save some 
space, put 

Xi:=L(r,_i,ri]-c(ri-ri_i). 

Donsker's theorem says that 



—^Y^Xun>^\ ^{Bu,u>0), 
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weakly in D[0,oo), as n — )■ oo, where S is a standard Brownian motion. Let 

ipn{t):=^, t>0. 
n 

Since (/?„ converges weakly, as n — ^ oo, to the deterministic function (At, i > 0) and since 
composition is a continuous operation, the continuous mapping theorem tells us that 

f — ^ Y. Xi.u>^\^ {Bxu,u > 0) = X^/^B, 

and this readily implies that the last term in (I30p converges weakly to a Brownian motion. 

D 

It is now easy to see how the quantity T[i,j], the maximum length of all paths from 
i to j, behaves. A sufficient condition for T[i,j] to be positive is that there is a skeleton 
point between i and j. Therefore, keeping i fixed, the probability that eventually for all j 
sufficiently large T[i,j] > is at least equal to the probability that eventually there is a 
skeleton point in [i,j], and this is certainly equal to one. So, eventually, any two points are 
connected, a.s. 

Moreover, 

r[r„r,] = L[r„r,]. 

Indeed, Fj is connected to every larger vertex and any vertex smaller than Tj is connected 
to Tj. Thus, if a path from some u > Fj to some v < Tj has length L[rj,rj] we necessarily 
have u = Ti and v = Tj and this shows the equality of the last display. 

If n is large enough so that there is at least one skeleton point in [0,n], we have that 
~^ n and 

L[ri,r$j<r[o,n]<L[ro,r$„+i], 

where ^n is the number of skeleton points in [0, n]. Therefore we immediately obtain: 
Theorem 3. // [CI] and [C2] hold then T[0,n\/n — )• C, as n ^ oo, a.s. 

Same rationale shows: 
Theorem 4. Suppose [CI] and [C3] hold. Then Theorenn\^ holds with T in place of L. 

8 Directed slab graph 

Recall that we started with vertex set V = 'L and introduced a random partial order ~-^ 
by means of a random directed graph: 

i -^ j lii < j and 3 i = i^ < ii < ■ ■ ■ < i^ = j such that aio,Ji = • • • = ai^.j j = 1. (31) 

A natural generalisation is to replace the total order < on the vertex set 1/ by a partial 
order -< and substitute the i < j requirement in (j3ip above by the requirement that i ~< j. 
We here provide an example of such a generalisation. A major role in our analysis has been 
played by the assumption that the underlying probability measure is invariant by some shift 
9. Our example will also satisfy this assumption. 
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Let (I, :<) be a finite partially ordered set. We assume that / has a minimum and a 
maximum, denoted by and M, respectively. In other words, for all i,j, k £ I, 

(a) O^i^i^ M, 

(6) if i ^ J ^ i then i = j, 

(c) ii i ^ j ^ k then i <k. 

Consider y = Z x /. We call this vertex set a cylinder. In the case / = {0, 1, . . . , M}, with 
the usual ordering, we call V a slab. Elements of V will be denoted by (x,i), {y,j), etc. We 
introduce the component- wise partial ordering ■< on V by 

{x,i) < (y,j) ^^ (x,i) / (y, j) and x < y, i ^ j, 

and write {y,j) >■ (x,i) for the same thing. Next, we assign an edge {{x,i), {y,j)) to each 
pair of vertices such that (x,i) <; (?/ij) with probability ry-x,i,j, independently from pair 
to pair. This is done by means of random variables atx,i),(y,j)'- 

IP(a(x,i),(s/j) = 0) = 1 - P(a(a.^i)^(yj) = -oo)ry_x^i^j. 

We shall make this more formal in the sequel. The problem is, again, the behaviour of a 
longest path from {x,i) to (y, j). This length is denoted by T[{x,i), (y, j)]. We also define 
L[(a;,z), (y, j)] to be the maximum length of all paths starting from some (x',i') >> (x,i) 
and ending at some {y\j') > {y-,j)- 

An appropriate probability space for the model is now described. Let 8 = {5x,ij,x £ 
1',i,j G /) be a collection of independent {— oo, l}-valued random variables with 

assuming that Vx^ij = if x < or if z ;^ j. Next, let S^^', x G Z be a collection of i.i.d. 
copies of 6. The probability space O. is defined to contain infinite vectors w = {S^^>,x S Z). 
In other words, i7 = ({— oo, i}^><^x^)2 with {— oo, ly^^^^^ be the space of values of each 
S^^' , and with P being a product measure. A shift on f] is taken to be the natural map 

uj = {x^d^''^)^euj = {x^S^^+^^). (32) 

Clearly, P is preserved by 6. The random variables a[x,i),{y,j) are now given by 

Oi{x,i),{y,j)[^) =Oy-x,i,j 

and it is easy to check their 6'-compatibility: a(^x,i),{y,j){(^'^) = Oi{x+i,i),(y+i,j)i^) ■ 
We introduce the following assumptions on the probabilities rx,ij- 

[DO] rx,i,i =: Px for sdli £ I 

[Dl] < pi < 1 

oo 

[D2] ^(l-pi)...(l-p,)<oo 

x=l 



[D2'] ^a;(l-pi)---(l-p,)<oo 



x=l 



[D3] For all i,j £ I with i -< j, we have ro^ij > 0. 
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Of these, the last one is not an essential condition. It is only introduced for convenience. 
We will comment on it later. Of course, [D2'] is stronger than [D2] and it will be used for 
the proof of the CLT. 

8.1 The random graph G[x,y\ 

The random directed graph G = {V, E) with V = "Z x I and E consisting of all 
((x,i), (y, j)) such that a(a;,i),(i/,j) = 1 is now a well-defined object. Let G[x,y] be the 
restriction of G on the vertex set [x, y] x I where x < y are two integers. Let 

L[x,y]:= rnax L[{x' , i) , {y' , j)] 

x<x <y <y 

be the maximum length of all paths in G[x,y]. We have ^-compatibility 

L[x,y]oe = L[x + l,y + l], 

and, by an argument analogous to the one used to obtain Q, we have the subadditivity 
property 

L[x, z] < L[x, y] + L[y, z]+l, x < y < z. 

Therefore, 

Ln/N := L[0,N]/N -^ G, as n -^ oo, a.s., 

for some deterministic constant C which, under the assumption [D2], is positive. 

8.2 The random graph G*^*) 

Let G^"^' be the restriction of G on the vertex set V x {i}, i ^ I. It is clear that each G^'^' 
is a line model as studied earlier. In fact, the G^^\ i ^ I are i.i.d. We denote by L^^>[x,y\ 
the maximum length of all paths of G^'*"' from some vertex x' > a; to some vertex y' < y. 
We shall let ^(*) be the skeleton of G^*). Then, assuming [Dl] and [D2], each ,5^(*) forms 
a stationary renewal process with nontrivial rate. Moreover, [Dl] implies that this renewal 
process is aperiodic. 

9 Central limit theorem for the directed cylinder graph 

We first describe the limiting process. To do this, we need the following. First, let 
{B^^\t),t > 0), i G /, be i.i.d. standard Brownian motions, all starting from 0. Second, let 
Hil, ^) be the Hasse diagram [16] corresponding to the partially ordered set /. This is a 
directed graph with vertex set I and an edge from i to j if there is no k, distinct from i 
and j, such that i ^ k ^ j. Let l = {lq, li, . . . ,Lr) be a path in H{I, <) starting from lq = 
and ending at tr = M. The length of the path is r = \l\. For each such path t, define the 
stochastic process (Z(t)i,t > 0) by: 

Z{t)t : sup {s('«)(to) + [B^'^Kh) - B^'^\to)] +■■■ + [i?(^l'')(tH) - i?(^l'i)(ti.|_i)]} 

0<to<ti<---<t|,|=i 

(33) 
and then let 

Zt:= max B{i)t, (34) 
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where the maximum is taken over all paths l from the minimum to the maximum element 
in the Hasse diagram. 

The main theorem of this section is as follows: 

Theorem 5. Let G he a directed cylinder graph and assume that [DO], [Dl], [D2'], [D3] 
hold. Let Ln be the maximum length of all paths in G[0,n]. There exists a constant k > 
such that 



init) ■■-- 



L, 



nt] 



Cnt 



t > 0, nGN 



K\ n 



converges weakly, as n —)■ oo, in the Skorokhod space L'[0,oo) equipped with the topology of 
uniform convergence on compacta, to the stochastic process Z defined in 



Proof. Since the =5^*-*^ i G / are independent aperiodic renewal processes, we have that 

o5^ := {re G Z : x E njg/=5^^*\ OL{x,i),{x,j) = 1 fo^' ^^ hj^l with i -< j} 

is also a renewal process. Indeed, Lindvall |26] shows that HiJ^^^^' is a stationary renewal 
process. Now o5^ is obtained from Hi^j^^^' by a further independent thinning with positive 
probability due to the convenient assumption [D3]. Condition [D2] implies that the rate of 
each =5^(*) is positive and this implies that the rate of njg/=5^(*^ is positive. Hence the rate 
of J^ is also positive. Call this rate A. We have < A < 1. Moreover, J/" is stationary: 
J/'oO = S^ . Enumerate now the points of 5^ by 



< r_i < To < < Ti < Ta < • • • 



We have E(r2 - Ti) = 1/A. If 



$„:=|^n[0,n]|, 



we have lim„^oo ^nln = A, a.s. Furthermore, C = AEL[ri,r2] < 1. Condition [D2'] implies 
that E(r2 - Ti)^ < oo and hence EL^*) [T2 - Ti]^ < 00. By Corollary [U the random variables 

(r2-ri, L«[ri,r2]), (r3-r2, L»[r3,r2]), ... 

are i.i.d., and since <y is obtained by independent thinning of nig/=5^(*', we further have 
that the rows of the last display are also independent when i ranges in /. 



[N,M] 




[0,0] '"i 



Figure 2: The skeleton for the slab graph and a longest path. 
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Consider next a path t = (to, ti, ■ ■ ■ , i-r), of length \i\ = r in the Hasse diagram H{I, :<) 
and define the quantities 

L*ii)r,:= max {L('^^[T,,Tj,] + L('-^^[T,„T,,] + ■ ■ ■ + lM[T,,._,,T,^]} 

l<JO<Ji<--<Jr=n 

L* := maxL*(i)„, 

I, 

where the last maximum is taken over all paths t from the minimum to the maximum 
element of the Hasse diagram. 

We now argue that the quantity of interest L„ is of order L* + Od(l) when n is large by 
providing an upper and a lower bound. The key observation is that when n is large, the 
number of points Tj < n grows at a positive rate (and hence to infinity). At each of these 
points, say Tj, the graph G[rj,rj] (being a vertical slice of G-see Figure 2) is precisely the 
Hasse diagram: 

G[rj,r,] = H{i,^), jez. 

Fix l' -< i" in /. Since Tj is a point in the skeleton of G'' ' , any x < Tj is connected to Tj in 
G^'' \ Similarly, Tj is connected to any y in G^'' K Since l' is connected to i" in G[rj,rj], it 
follows that, almost surely, there is path in G from any (x, l') to any (y, l"), ii x <Tj < y 
for some Tj G =5^ and if l' -< l". 

Assume that <^„ > 2. Let l = {lq, li, . . . ,Lr) be a path in H(I, :<) with lq = 0, Lr = M 
and consider integers 

1 < Jo < il < • • • < ir-l < jr = ^n- (35) 

Keep in mind that 

r$„ < n. 

By the construction of the set J^, the following is true: 

(ri,o) = (ri,io) -^ (rjo,io) -^ (rj,pii) -^ (rii,ii) -->••• -w {Tj^^^,ir) -> {Tj^,ir) = (r$„,M), 

where {x, l') ~^ (y, l") means that there is a path from (x, l') to (y, t") in G. Therefore 

Ln > L^'^'^TuT,,] + L(^i)[r,„,r,j + • • • + L(^'-)[r,,_„r,v], 

because the right-hand side is a lower bound on the length of the specific path chosen in the 
last display. By keeping t fixed and maximising over the jo , . . . , jr satisfying psp we obtain 
Ln ^ L*{i)n, and by maximising over l we obtain the lower bound 

Ln ^ -^<j>„- 

To obtain an upper bound, let vr* be a path that achieves the maximum in L.„. Assume 
that $„ > 1 so that, by the key observation above, (0,0) is connected to {n,M) in G. See 
Figure 3. Hence vr* is necessarily a path from (0,0) to (n,M). 

Let 

= Lq ^ ii ^ ■ ■ ■ ^ is = M 

be the distinct values of the /-components of the eleinents of vr* in order of appearance 
in TT*. (The sequence {lo,li, . . . ,Ls) is not necessarily a path in H{I,^).) So for each 
k = 0, . . . ,s — 1, there are vertices {xk, ik), ivk, ^fc+i) which are consecutive in the path vr*. 
Hence 
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Figure 3: Construction used in obtaining the upper hound. 



Xk<yk< Xk+i, for all /c = 0, 1, . . . , s - 1, 

where, by convention, we set Xs = n. The point of S^ prior to x^ is r$^ and, since vr* has 
maximum length, (r$^ , t^) is an element of vr*. By the maximality of vr* again, we have 
that Xk and yt are contained between two successive points of y (otherwise we would be 
able to strictly increase the length of the path) . Hence 



r$, <Xk< Vk < ri+<j.,. < 2:^+1, for all /c = 0, 1, . . . , s - 1. 



(36) 



We thus have 



K = \n*\=L^^^^[0,T^]+L^^^^[T,,T^J+^U^^''\T^^^,Xk]+l+L^'^+^Hyk,r^^^^ 

fc=o "^ ^ 



n\ 



Due to (I36l) . we have 
Moreover, 

L('"^[o,ri]<L('°)[ro,ri] 



(37) 
(38) 

(39) 
(40) 



Each of the right-hand sides of (f37|) . p9|) and (|iO]) is bounded above by maxo<j<<j>„ L^') [Fj, Fi+j] 
If we then define 

and use (l38l). we obtain 



i^n < Cn + M + 5^ l(^'=+i) [r$,^ , r$,^^ J . 



fe=0 



Since for each sequence = io ^ '-i ^ " " " ^ '-« = -^ of distinct ordered elements of / we can 
find a path in the Hasse diagram containing these elements, it follows easily that 

Ln<Cn + M + Ll^, 
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which gives the upper bound. The upper bound is close to L^ in the sense that the sequence 
the Cn are of order 1 in distribution, i.e. that ((^„) is tight random sequence. On the other 
hand, nt = r$, j, — Fi + Od(l). It is thus clear that the weak limit of in and that of 



K\/n 
if it exists, will be identical. Setting 



«n [t) ■= 7= ' V5n(ij :- 



we have 

^:(i)=C(^n(t)), (41) 

and so the weak limit of ^* is equal to that of ^J^ (if this exists) composed by the function 
{Xt}. 

To show that the weak limit of £* exists and find it, define the function ip : -D[0, oo) — )■ 
Z)[0,oo)by 

V'(/3«,iG/)(t) :max sup |/3('"^(to) + [/3^'^^(ii) - /3^''^(to)] + • • • 

^ 0<to<<l<-<<r=t I 

\L\=r 

•••+[/3(^'-)(i,)-/3(^'-)(t,_i)] 

where the maximum is taken over all paths l from the minimum to the maximum element 
in the Hasse diagram H{I, :<). The function ip is continuous (with respect to the topology 
of uniform convergence) . Let 

4-^(t):= '^^'^^'^^''^'"^'^"f^^'^'"^'"'^^\ t>0, ZG/, 

where 

a2:=var{L«[ri,r2]-C7(r2-ri)} 

Since L'^*'[rj,rj+i], j > I, i ^ I are i.i.d. with common variance u^ we have (Theorem [2]) 
that 

(s«,iG/)^(i3»,iG/) (42) 

where i?*^*', i G / are i.i.d. standard Brownian motions. Let 

K := A^V 

and observe that 

C(t) = A-i/2.V^(.«,zG/)(t). 

By ()42p and the invariance principle. 

By the relation (I4ip and the remark following it, we have 

and the right-hand side is equal in distribution to Z (defined by (f33l) - (f3ll) ). D 

The remarks at the end of Section [7] also apply in the current case. We can easily conclude 
that Tn, the maximum length of all paths from (0,0) to {n,M), has the same asymptotics 
as L„. In particular. Theorem [5] holds if we replace L„ by T„. 
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10 Connection to last passage percolation 

Consider now the case 

/ = {0,1,...,M} 

with the usual ordering. Assumption [D3] can be substituted by 

[D3'] For all 1 < i < M we have ro,i_i,i > 0. 

Let Gm be the corresponding random directed cylinder graph, referred to as slab graph here. 
In particular, we can think of Gm as the restriction of a graph Goo on the vertex set Z x Z4. 
where two vertices (x,i) and {y,j), with {x,i) <^ {y,j), are connected with probability 
Py-xj-i that depends on the relative position of the two vertices on the 2-dimensional 
lattice. 

The problem here becomes that of a last passage percolation , although the model is not 
the standard nearest-neighbour one. Physically, we can think of tunnels which run upwards 
(or in directions southwest to northeast) and fluid moving in tunnels. It takes one unit of 
time to cross a specific tunnel. We are interested in the particle that starts from (0, 0) and 
reaches (n, M) in the largest possible time. Since the Hasse diagram of the set {0, 1, . . . , M} 
with the natural ordering is the linear graph with edges from i — 1 to i, !<«< M, the 
limit process Z is given by the simplified expression 

Zi = _ max {B^''\to) + [B'^'\h)-B'^'\to)] + --- + [B'^^'HtM)-B^^'HtM-i)]}, t>0. 

The latter process is a Brownian last passage percolation process. As was shown in [6ll2H[32] 
it is a non-Gaussian process with marginal distribution 

Zt = Vt ■ Am, 

for each t > 0, where Am is the largest eigenvalue of a random (M + 1) x (M + 1) Gaussian 
Unitary Ensemble (GUE) [29] . 

Tracy and Widom |341 I35j showed that, as M — )• 00, the following weak limit holds: 

m1/6(Am - 2\/M) => Ftw, 

with Ftw being the Tracy- Widom distribution whose hazard rate equals f^ q{x)'^dx, where 
q{x) satisfies a Painleve II equation; see [21 eq. (3.1.7)]. For an account on the universality 
of this distribution, see, e.g., [TTj. A number of interesting results have been proved relating 
this limiting distribution with certain stochastic models. These models include longest 
increasing subsequence |5], last passage percolation, non-colliding particles, tandem queues 
[HI [21], and random tilings [24] ■ For the last passage percolation, in particular, this limit 
is known to appear in two cases. The first is the Brownian last passage percolation. The 
second is the last passage percolation model with exponential (or geometric) weights. In 
this model one puts independent and identically distributed exponential random variables 
in the vertices of Z^ and considers the maximum L{M, N) of the sums of the weights over 
all directed paths from (0, 0) to (M, N). It was shown in [23] the random variable L{N, N), 
properly normalized, converges to the Tracy- Widom distribution as N goes to infinity. In 
[8], more general weights were considered and an analogous result for the random variable 
L{N,N"') (for an appropriate a depending on moment conditions) was obtained. It is then 
natural to conjecture that a similar phenomenon occurs in our slab graph too. 
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